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■q^-, Abstract 

O ! This study revisits the problem of advective transfer and spectra of a diffusive scalar field 

^ ! in large-scale incompressible flows in the presence of a (large-scale) source. By "large-scale" it 
{3jq! is meant that the spectral support of the flows is confined to the wave-number region k < k^, 
5 ! where kd is relatively small compared with the diffusion wave number k^. Such flows mediate 
couplings between neighbouring wave numbers within k^ of each other only. It is found that 
the spectral rate of transport (flux) of scalar variance across a high wave number k > kd is 
bounded from above by UkdkQ{k,t), where U denotes the maximum fluid velocity and Q{k,t) 
Ch ! is the spectrum of the scalar variance, defined as its average over the shell {k — kd,k + kd). For 
^ ! a given flux, say ^ > 0, across k > kd, this bound requires 

^: Q{k,t)>-^k-\ 

^. Ukd 

o , _^ 

, This is consistent with recent numerical studies and with Batchelor's theory that predicts a k 

^ 1 spectrum (with a slightly different proportionality constant) for the viscous-convective range, 

P^i which could be identified with {kd, k,^). Thus, Batchelor's formula for the variance spectrum is 
recovered by the present method in the form of a critical lower bound. The present result applies 

, to a broad range of large-scale advection problems in space dimensions > 2, including some filter 

O ! models of turbulence, for which the turbulent velocity field is advected by a smoothed version 

I of itself. For this case, Q{k,t) and "i? are the kinetic energy spectrum and flux, respectively. 
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1 Introduction 

The problem of scalar transport and mixing in turbulent fluid flows has been a subject of active 
research for decades, dating back to the late 1940s. Early studies by Obukhov [1] and Corrsin [2] 
applied Kolmogorov's theory of turbulence in a straightforward manner. They found that the 
scalar (fluid temperature in their case) variance behaved in the same manner as the turbulent 
kinetic energy, cascading via a k~^^^ range to a diffusion range at high wave numbers k for 
disposal. This result is supposed to apply to cases of relatively small diffusivity k and viscosity 
v in the regime k ^ u, for which the viscous dissipation and diffusion ranges coincide. Batchelor 
[3] considered turbulent flows at moderate Reynolds numbers in the regime of large Prandtl 
or Schmidt number Pr = u/k ^ 1, for which there exists a broad viscous-convective range 
ky <^ k <^ kfi between the viscous dissipation wave number k,j and diffusion wave number fc^- 
He found that in this range, the scalar variance spectrum F{k) scales as k~^ and is given by 

F{k) = ^k-\ (1) 
7 



where x is the rate at which the scalar variance is dissipated, i.e. the spectral rate of variance 
transport or variance flux, and 7 is an effective least-ratc-of-strain parameter given by 7 = 
C(e/z/)^/^. Here e denotes the mean rate of kinetic energy dissipation and C is a constant of 
order unity. From the Obukhov-Corrsin and Batchelor theories one may visualize a picture of 
scalar advection in flows at moderate Reynolds numbers in the limit of large Pr, in which a 
hybrid spectrum obeys the Obukhov-Corrsin k~^^^ scahng in the fluid inertial range followed 
by the Batchelor scahng in the viscous-convective range [4]. These pioneering theories have 
been considered to be breakthroughs and attracted considerable interest to the subject during 
its infancy [5-9]. Recently, fundamental issues in geophysical, environmental and industrial 
applications have sparked a surge in the area, resulting in a huge body of research [10-30] on 
a variety of dynamical aspects. Another reason for this surge is that computers have become 
increasingly capable of taking on a scientific problem of this magnitude. Within the past few 
years, numerical evidence in support of the Batchelor theory and its predicted spectrum 
has accumulated considerably [4, 22, 30]. However, this is far from conclusive as the viscous- 
convective ranges accessible to modern computers are still quite limited. Furthermore, a number 
of studies [16, 24-26] have either argued for or found spectra shallower than the Batchelor 
spectrum. For these reasons, as well as the phenomenological nature of the Obukhov-Corrsin 
and Batchelor theories, further theoretical consideration and numerical analysis (with higher 
resolutions whenever possible) continue to be desirable. 

In this study, we revisit the advection-diffusion problem, by carrying out a simple but 
rigorous analysis of the advective transfer term leading to a conclusion that is consistent with the 
Batchelor picture [3] and with recent numerical results [4, 22, 30]. We consider large-scale flows, 
meaning that the tail of the Fourier representation of the flows beyond some flnite wave number 
kd is identically zero or at least can be ignored. Such smooth flows are relevant for practical 
purposes as most advection-diffusion problems in the geophysical and environmental contexts 
are primarily concerned with large-scale advecting flows. They may even model Navier-Stokes 
turbulence at moderate Reynolds numbers if kd belongs to the viscous dissipation range and if 
the exponentially decaying tail of the velocity flelds beyond kd can be ignored. These large-scale 
flows can mediate transfer of the scalar variance between neighbouring wave numbers within 
kd of each other only. It is found that the variance flux across a high wave number k > kd is 
bounded from above by UkdkQ{k, t), where U denotes the maximum fluid velocity and Q{k, t) is 
the spectrum of the scalar variance, deflned as its average over the shell {k — kd,k + kd). From 
this upper bound, it can be readily deduced that in the high k limit, the flux diminishes if 
Q{k, t) becomes steeper than k~^. Hence, no spectra steeper than k"^ could support a nonzero 
variance flux to the vicinity of the diffusion wave number A;^ in the limit of large ki^/kd, which 
may be identifled with Pr (or P^^"^) [31]. Given a persistent source and in the limit k^/kd 00, 
this result implies a divergence of the scalar variance not slower than logarithmic in k even if a 
variance cascade to the small scales is realizable. The present results apply to both passive and 
active scalars in large-scale flows, provided that in the active case, the excitation of the wave 
numbers k > kd oi the flows by nonlinear feedback mechanism can be ignored. They also apply 
to fllter models of turbulence, for which the full turbulent velocity is advected by a smoothed 
version of itself. In this case, Q{k, t) is replaced by the kinetic energy spectrum, and the flux 
in question is the kinetic energy flux. 



2 Preliminaries 



In this section, we briefly describe the advection-diffusion equations in spectral form, principally 
to illustrate the coupling locality, which plays a key role in the present analysis. We then recall 
the variance conservation law and set out a few notations employed in this paper. 

The advection-diffusion equations governing the evolution of a diffusive field 9{x, t) advected 
by incompressible flows u{x,t) are 

BO 

— + u-ve = KAe + f, V • w = 0, (2) 

where k, is the diffusivity and f{x,t) is a (large-scale) source. The spectral support of u{x, t) is 
assumed to be confined to the region k < kd, where ka is a finite wave number. We consider Eq. 
(2) in an n-dimensional (n > 2) periodic domain, enabling us to express our results conveniently 
in terms of spatial averages of dynamical quantities. These results can be seen to carry over 
to an unbounded domain with minimal change. All fields are assumed to have zero spatial 
average. The advected field 9{x,t) can be either passive or active. In the latter case, the 
nonlinear feedback mechanism by 9{x,t) on u{x,t) can be arbitrary, as long as it does not 
"irregularize" u{x,t) by exciting the small scales of u{x,t) corresponding to k > ka to the 
extent that these scales can no longer be ignored. Furthermore, 9{x,t) can be a vector, such 
as the fluid velocity in some filter models of turbulence (for which a pressure term is included) . 
The Fourier representations of u{x,t) and 9{x,t) are 

u{x,t) = ^ ■u(A;, t) expjifc • a;} (3) 

k=\k\<kd 

and 

9{x,t) = Y,^(k,t)e^^{ik-x}, (4) 

k 

respectively. Here fe 7^ is the wave vector and w(fc, t) and 9{k, t) are the Fourier transforms of 
u{x,t) and 9{x,t), respectively. The reality of u{x,t) and 9{x,t) requires u{k,t) — u*{—k,t) 
and 9{k,t) — 9*{—k,t). The incompressibility of u{x,t) further requires k ■ u{k,t) — 0. In 
spectral form, the first equation of Eq. (2) becomes 

^9{k,t)= Yl k' ■u{k",t)9{k',t)-Kk^9{k,t) + f{k,t), (5) 

where f{k,t) is the Fourier transform of f{x,t). The incompressibility of u{x,t) manifests 
itself in Eq. (5) through the fact that k' ■ u{k" ,t)9{k' ,t) = if fc' and k" are coUinear. The 
triad relation k — k' + k", together with the constraint k" = \k"\ < ka, imphes that k' — \k'\ 
satisfies \k — k'\ < kd- This means that a given k > ka can couple with other wave numbers 
within the shell {k — kd,k + kd) only. This coupling locality has a significant consequence as 
will be seen in the next section. 

Given a "reference" wave number A;* (A;* ^ kd), let us denote by 6'<, 6'>, 9i, and 9^ the 
components of 9 spectrally supported by the ball b = {k : k < k^}, its complement B = {k : k > 
k^}, the inner shell Si = {k : k^ — kd < k < A;*}, and the outer shell = {k : k^, < k < k^:-\-kd}i 
respectively (see Fig. 1). For example, 

9> {x, t)=J2 t) exp{ifc • x}. (6) 

fce-B 



For k > kd, let Q{k, t) denote the averaged variance taken over the shell {k — kd,k + ka). For 
example, Q{k*,t) is given by 



e(fc.*) = ^^^^, (T) 

where (■) denotes a spatial average. It is evident that Q{k,t) approximates the usual spectrum 
F{k,t). In general, the approximation can become increasingly better for higher k as the shell 
{k — kd,k + kd) becomes thinner, in the sense that the ratio of the shell thickness 2kd to its 
radius k becomes smaller. For power-law spectra, Q{k,t) actually approaches F{k,t) in the 
limit k/kd ^ oo. For example, for the Batchelor spectrum given by Eq. (1), we have 

2A;d7 Jk.-kd 2kd^ - kd 2kd^ \ k^-kdj 

which indeed tends to xil^*)^^ in the limit k^/kd oo. Thus Q{k,t) tends to F{k,t) in the 
limit of high k (at least for power-law spectra). We will express our results in terms of Q{k,t) 
instead of F{k,t) since the former arises more naturally in the present context. 

Given periodic functions (p and ip having zero mean and bounded mean-square gradients, 
we have 

{(f)u-Vilj) ^ -{i/ju-Vcj)). (9) 

This identity gives rise to a wealth of conservation laws, particularly the variance conservation 
law, and is used repeatedly in this study. 



3 Main results 

We now present the main results of this study. First, we elaborate on the locality of the 
variance transfer and then derive the lower bound for Q{k,t) as described earlier. Second, we 
show that given a persistent source, {6'^) grows without bound in the limit A;^ oo, irrespective 
of the underpinning dynamics. Third, we compare the present finding with a recent result [31] 
derived on the assumption of bounded velocity gradients, i.e. |Vw| < oo, which is a weaker 
condition than the present one. Finally, the slight discrepancy between the present finding and 
Batchelor's formula (1) is discussed. 



3.1 Bounds for variance flux and spectrum 

The governing equation for the evolution of the small-scale variance (6'^) is obtained by multi- 
plying Eq. (2) by 9y and taking the spatial average of the resulting equation, 

^|(^>) + '^{M>\') = -(^>« • v^) = -{e^u ■ v^<), (10) 

where Eq. (9) and the linearity of the advection term have been used and the forcing term 
vanishes as B is assumed to be source free. The triple-product term (fiux term) in Eq. (10) 
represents the net variance transfer across /c* into the region k > k^, which drives the small- 
scale dynamics. At the modal level, this fiux term consists of triple-product terms of the form 
9{k,t)e{k',t)k' ■ u{k\t), where fc e S, fc' e 6 and fe = fc' + k" . Since k" < kd, this triad 




Figure 1. A schematic description of the spectral supports for u and for the components 
9y, 9i, and 9e of 9. The flux term involves only the wave numbers within the shell (/c* — 
kd, + kd), which supports 9i and 9e- 



relation implies that k and k' can couple only if A; — A;' < kd- Hence, only modes in 9i and 9e, 
i.e. within the wave number shell (A;* — kd,k^ + kd), contribute to the flux term (see Fig. 1). 
For this reason, the variance transfer can be considered as being highly local, particularly at 
high k, where the shell {k^, — kd, /c* + kd) becomes relatively thin (radius becoming larger but 
thickness remaining fixed). Here, we use the term "highly local" to emphasize the fact that 
k/k' ^ 1 in the limit k*/kd —>■ oo. This term is to distinguish the present couplings from those 
of a lesser degree of locality between k ^ k' ^ k^ via k" ^ k^, where the ratio k/k' remains 
strictly greater than unity in the same limit. Such couplings are clearly absent from the flux 
term. From this geometric consideration, we can write 

(^>'u- V^<) = {9,u-V9i). (11) 

Substituting this result into Eq. (10) yields 

< uk^mY'^mY^^ < ^ {m') + (^5) = ukdk,Q{k,,t), 

where, as wc recall, U denotes the maximum fluid velocity and Q{k,t) is the variance spec- 
trum defined by Eq. (7). In Eq. (12), the Cauchy-Schwarz inequality and the self-explanatory 
(Poincare-type) inequality (|Vft,:p) < kl{9f) have been used. The bound for the flux term in 
Eq. (12) is interesting and can be readily interpreted in what follows. 

For a positive flux through k^,, say t?*, the final estimate in Eq. (12) implies that 

■J?* < Ukdk*Q{k*, t), 

or equivalently 

eik.,t)>4irK'. (14) 
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It follows that a positive /c-independent flux is possible only if Q{k,t) becomes no steeper than 
k^^ (pointwisc) for high k. This constraint is consistent with Batchelor's theory that predicts 
a k~^ spectrum for the viscous-convcctivc range, which could be identified with (k^, fc^)- Since 
Eq. (14) implies a divergence of (9^) toward the small scales at least as rapid as logarithmic 
in k, a positive variance flux to ever smaller scales (including those that diminish no more 
rapidly than (Ink)^^) requires a priori an unbounded variance "passage". This is in a sharp 
contrast to the classical direct energy cascade (and the Obukhov-Corrsin variance cascade), 
which is supposed to proceed through an inertial range virtually free of energy. In some sense, 
the energy cascade is rather "rushing" , whereas the variance cascade of the present case (if 
realizable) would be far less dramatic, "leaking" through a fully fllled inertial range. In the 
presence of a persistent scalar source, (6'^) necessarily grows without bound in the limit k^ ^ oo. 
for obvious reasons. On the one hand, a variance cascade to ever-smaller scales already requires 
at least a logarithmic divergence of {6"^) toward the small scales. On the other hand, if such a 
cascade is unrealizable, the injected variance is necessarily trapped at the large scales, thereby 
resulting in their unbounded growth. For the sake of completeness, this argument will be made 
more quantitative in the next subsection. 

3.2 Unbounded variance growth in the limit ^ oo 

Similar to Eq. (10), the governing equation for the evolution of the large-scale variance (^|) is 
obtained by multiplying Eq. (2) by 6'< and taking the spatial average of the resulting equation, 

= {0eU ■ VOi) + 1? > -Ukdk,e{K, t) + (15) 

where Eqs. (9) and (11) have been used and the inequality is a straightforward application of 
the upper bound for the flux term derived in Eq. (12). In Eq. (15), ■»? = (^</) is the scalar 
variance injection rate. For some large time t — T, say T — 1/(2«;A;^), let Q denote the average 
over [0,T] of a dynamical quantity Q. Taking the time average of Eq. (15) and re-arranging 
the terms in the resulting equation yields 

nkl{el) + n{\ve<\^) > ^ - ukdk,e{k,), (i6) 

where the initial value of {9^) has been omitted for convenience. Upon making the substitution 
(|V^<P) < kl{el) in Eq. (16), we obtain 

nkl{9l) + nkl'{9^>d-UkdkMk*). (17) 

In accord with a persistent source, let us assume i) > 0. Now in the limit k^ ^ oo (k ^ 0), if 
there exists no /c* < oo such that the right-hand side of Eq. (17) is positive, then (9^) diverges 
toward the small scales as discussed above. On the other hand, if there exists k^ < oo such that 
the right-hand side of Eq. (17) is positive, then ((^|) + (^<)) — > oo. It follows that (^|) — > oo, 
and hence (9^) — > oo. Thus, {9'^) diverges regardless of whether or not there is a variance 
cascade. 

3.3 Discussion 

When u{x,t) is not restricted to the large scales, there are no constraints on k' and k in the 
triple-product terms 9{k, t)9{k' , t)k' ■ u{k", t) contributing to {9>u • V^<). The flux term then 



involves, in principle, conplings for every k' < and k > /c*. The presence of nonlocal couplings 
(between k' ^ k^ and k ^ A;* via k" ^ k) and the other type of local couplings (between 
k ^ k/ ^ k^: via k" ^ k^,) mentioned earlier effectively makes the flux term unmanageable 
by the present method, in the sense that its analytic estimates would be too excessive for 
meaningful interpretations. For this case, Tran [31] finds by a different approach that if the 
advecting velocity fields have bounded gradients, then diffusion anomaly, i.e. a variance cascade 
to ever smaller scales, requires the variance at the small scales to be no less than that provided 
by the Batchelor k"^ spectrum. This constraint is weaker than the present one as it does not 
rule out the possibility of bounded variance corresponding to non-power-law spectra having 
gaps of severe variance deficiency in the intermediate wave-number region, provided that the 
variance requirement at the small scales is met. The present finding, by exploiting the high 
locality of the variance transfer for large-scale advecting flows, rules out this possibility. The 
variance is required to grow without bound either via bounded spectra not steeper than 
(pointwise) if a variance cascade is realizable or via unbounded spectra if otherwise. 

In the absence of a scalar source, a finite variance reservoir cannot support a A;-independent 
flux because such a flux requires an unbounded variance "passage" as we have concluded. Our 
result allows for no significant "chunk" of a given initial variance reservoir (^q) < oo at large 
scales to break away and cascade to the small scales by itself. Rather, it suggests a gradual 
spread out of (^^q) ever more thinly in wave-number space, giving rise to a diminishing flux, 
which can be readily estimated. Suppose that at a later time, a k~^ range gets established 
from kd to A;* ^ ka or beyond. Then, in this range, the spectrum Q{k,t) is bounded by 
Q{k,t) < {9l)k~^/ \n{k^/kd). Upon substituting this into Eq. (13), we obtain 



This means that diminishes at least as rapidly as {\n{k^/ka))~^. Note that although a 
logarithmic decay of the flux can be expected on heuristic grounds, Eq. (18) may not be 
rigorously derived without the constraint (13). 

The present bound (14) for Q{k,t) resembles the Batchelor formula (1) in every aspect 
except that Ukd in Eq. (14) plays the role of 7 in Eq. (1). This apparent discrepancy, however, 
can be reconciled if we reformulate the present problem in accord with the Batchelor setting. It 
can be seen that the product Ukd is essentially an upper bound for the velocity gradients |Vii|. 
So if we identify v{Ukci)'^ ~ z/|Vttp with the kinetic energy dissipation rate e in the Batchelor 
setting of turbulent advection, then we obtain Uk^ ~ (e/z/)^/^ ^ 7. Hence, Eqs. (1) and (14) 
agree. This is no surprise because the Batchelor problem would reduce to the present case upon 
the hypothesis that the exponentially decaying tail (beyond ka) of the turbulent velocity field 
contributes negligibly to the advective variance transfer. 

As it stands, Eq. (14) captures the intuitive physical fact that for fixed U , flows at larger 
scales (smaller kd) are poorer transporters as scalar spectra having larger spectral amplitudes, 
i.e. larger factors ^*/{Ukd)i would be required to support the variance flux across /c*. 



4 Turtulent energy transfer by large-scale advection 




(18) 



h\{K/kd) ' 



The above results apply to the energy transfer by large-scale advection in turbulence. Namely, 
the advection of the turbulent velocity by its large-scale component alone results in a contribut- 



ing energy flux that vanishes at high k if the energy spectrum becomes steeper than k~^. For 
the Kolmogorov k^^^^ spectrum, this means that the large-scale advection contributes negligi- 
bly to the direct energy transfer. On physical grounds, this is consistent with the expectation 
that the large scales, while advecting the turbulent eddies, do not stretch them significantly. 
The remaining of this paper is devoted to detailed elaboration of this fact. 
We begin by recalling the Navier-Stokes equations 

dv 

— + (v • V)v + Vp = i/Av + /, V • -y = 0, (19) 

where v{x^t) is the fluid velocity, p{x,t) is the pressure and f{x,t) is a large-scale forcing. Let 
w be a large-scale component of v, as defined by Eq. (3), and u' be its small-scale complement, 
i.e. V — u + u'. Furthermore, let v^, Vy, Vi, Vg and V{k) be defined in the same ways as 
^i: ^ud ©(/c), respectively. Note that V{k) is approximately twice the usual kinetic 
energy spectrum and that all the components of v so defined are incompressible. Similar to 
Eq. (12), the evolution of the small-scale energy (|'u>p)/2 is governed by 

^^(l^>r) + i^(|Vt^>r) = -{v> ■ {v ■ V)v<) = -{v> ■ {u ■ V)t;<) - (v> • {u' ■ V)v<) 

= -{ve ■ {u ■ V)Vi) - {v> ■ {u' ■ V)v<) 

< UkdKV{K, t) - (v> • (w' • V)v<), (20) 

where the forcing and pressure terms vanish as the region under consideration is assumed to 
be force free and v> is incompressible. In the final equation of Eq. (20), the first term on the 
right-hand side represents an upper bound for the energy transfer across /c* due to large-scale 
advection and the second term is the energy transfer across k^ due to small-scale advection. 
The former vanishes for high k^ if V{k,t) becomes steeper than k~^. This means that the 
latter is solely responsible for the direct energy cascade in the classical picture of turbulence, 
for which the k~^^^ energy inertial range is far too steep for the former to make a non-negligible 
contribution. At the modal level, this result is consistent with the expectation that triad 
interactions involving well-separated scales (those due to large-scale advection in the flux term 
{Ve-{u-'V)vi)) are relatively weak. Note that not all triads of well-separated scales are contained 
within ■ {u ■ V)t;j) as the term (v> • {u' ■ V)t;<) also has this type of triads. Such triads are 
formed by large-scale modes in v< and small-scale modes in v> and u'. Like their counterparts 
in {ve ■ {u ■ V)vj), these can be shown to be relatively weak and not responsible for the classical 
direct energy cascade. 

The above result may be applicable to models of turbulence that are derived by regularizing 
the Navier-Stokcs equations by a variety of filtering techniques (sec Graham et al. [32] for 
a discussion of several such models). For example, let us consider the Lcray model, obtained 
by dropping the small-scale component u' from the advecting velocity in the Navier-Stokes 
system, i.e. 

dv 

— + {u- V)v + Vp = i/Av + /, V • v = 0. (21) 

at 

For this simple model, the governing equation for (|i>>p)/2 is given by Eq. (20) without the 
small-scale advection term (t;> • {u' ■ V)v<). As a consequence, the classical direct energy 
cascade is not realizable for the reason discussed in the preceding paragraph. Instead, the 
energy behaves in the same manner as the variance (9^) described earlier. Namely, the energy 



either cascades to the small scales via spectra not steeper than or else accumulates at 
the large scales. Given a persistent source of energy, i.e. {v ■ f) > 0, the energy necessarily 
grows without bound in the inviscid limit. Equation (21) resembles regularization models of 
turbulence, which have been studied widely as alternatives to subgrid-scale models [32] and for 
which the k^^ scaling for the energy spectrum has been found by phenomenological arguments. 
The present result provides a different perspective to this possible scahng. 

In passing, we would like to note that the question of realizability of a (variance or energy) 
cascade and the associated k"^ (or shallower) spectrum cannot be resolved by the present 
analysis. This question is challenging because a lower bound for the flux term is highly infeasible, 
even for very simple flows. Given this difficulty, one may be better off resorting to numerical 
methods. What we have shown here is that if there is a cascade, then it must proceed through 
spectra not steeper than k^^ (pointwise). The critical k^^ scaling can be seen as most plausible 
for a number of reasons. In particular, it would correspond to a cascade of maximal spectral 
extent. 

5 Concluding remarks 

In summary, we have examined the advective transfer and spectral scaling of a diffusive field 
9{x,t) in large-scale incompressible flows u{x,t), whose spectral support is conflned to the 
wave-number region k < kd, for some finite wave number kd, which is relatively small as 
compared with the diffusion wavenumber k^.- The main result obtained is the upper bound 
Ukdk^,Q{k^,,t) for the variance fiux across a high wave number k^ > k^. Here U denotes the 
maximum fiuid velocity and Q{k,t) is the variance spectrum, defined as its average over the 
shell {k — kd, k + kd). The derivation of this bound exploits the very fact that the advecting fiows 
under consideration mediate variance transfer between neighbouring wave numbers within kd 
of each other only. The derived bound implies that for k ':$> kd, Si nonzero /c-independent flux 
is possible only if Q{k,t) becomes no steeper than k~^ (pointwise). This result is consistent 
with Batchelor's theory and with recent numerical and theoretical results [4, 22, 30, 31]. A 
novel element of the present findings is the pointwise constraint on Q(k,t) in Eq. (14). Given 
this constraint and a persistent source, the variance is required to grow without bound in the 
limit ki^/kd oo {k, 0), irrespective of the underpinning dynamics. The present results 
have been shown to apply to the Leray model of turbulence, for which the turbulent velocity 
is advected by a smoothed version of itself. Furthermore, they apply equally to different space 
dimensions > 2 as the analysis is dimension- independent. Finally, the possible relation between 
the advecting and advected fields does not enter the calculations. Hence, the results are vafid 
for both passive and active scalars, provided that in the latter case, the nonlinear feedback 
mechanisms on the fiows do not alter their large-scale designation. 
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